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Abstract 

We review the results obtained in an Effective Hamiltonian (EH) approacli for the 
three-quark systems. The EH is derived starting from the Feynman-Schwinger 
representation for the gauge-invariant Green function of the three quarks propa- 
gating in the nonperturbative QCD vacuum and assuming the minimal area law 
for the asymptotic of the Wilson loop. It furnishes the QCD consistent framework 
within which to study baryons. The EH has the form of the nonrelativistic three- 
quark Hamiltonian with the perturbative Coulomb-like and nonperturbative string 
interactions and the specific mass term. After outlining the approach, methods of 
calculations of the baryon eigenenergies and some simple applications are explained 
in details. With only two parameters: the string tension a = 0.15 GeV^ and the 
strong coupling constant = 0.39 a unified quantitative description of the ground 
state light and heavy baryons is achieved. The prediction of masses of the doubly 
heavy baryons not discovered yet are also given. In particular, a mass of 3660 MeV 
for the lightest baryon is found by employing the hyperspherical formalism to 
the three quark confining potential with the string junction. 
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1 Introduction 



The heavy flavor physics started with the discovery of the charm quark in 1974. It was 
followed by the discovery of the beauty quark in 1977 and of the top quark in 1995. Soon 
after the discovery of the charm quark several charmed baryons and mesons have been 
identified. The discovery of single charmed or beauty baryons has been followed by the 
discovery of the Be meson P] and recently by the extensive search of doubly charmed 
baryons [2]. 

Doubly heavy baryons are baryons that contain two heavy quarks, either cc, be, or bb. 
Their existence is a natural consequence of the quark model of hadrons, and it would 
be surprising if they did not exist. In particular, data from the BaBar and Belle col- 
laborations at the SLAC and KEK B-factories would be good places to look for doubly 
charmed baryons. Recently the SELEX, the charm hadropro duct ion experiment at Fer- 
milab, reported a narrow state at 3519 ± 1 MeV decaying in A'^K~n~^, consistent with 
the weak decay of the doubly charmed baryon [2^. The candidate is 6.3a signal. 

The SELEX result was recently critically discussed in Whether or not the state that 
SELEX reports turns out to be the first observation of doubly charmed baryons, studying 
their properties is important for a full understanding of the strong interaction between 
quarks. 

Estimations for the masses and spectra of the baryons containing two or more heavy 
quarks have been considered by many authors jlj. The purpose of this paper is to 
present a consistent treatment of the results of the calculation^ of the masses and wave 
functions of the doubly heavy baryons obtained in a simple approximation within the 
nonperturbative QCD. As a by-product, we also report the masses and wave functions 
for light and heavy baryons. 

The paper is organized as follows. In Section 2 we briefly review the Effective Hamiltonian 
method. In Section 3 we discuss the hyperspherical approach which is a very effective 
numerical tool to solve this Hamiltonian. In Section 4 our predictions for the ground-state 
spectra of light and heavy baryons are reported and a detailed comparison of spectra of 
doubly heavy baryons with the results of other approaches is given. Section 5 contains 
our conclusions. 



2 The Effective Hamiltonian in QCD 

Starting from the QCD Lagrangian and assuming the minimal area law for the asymptotic 
of the Wilson loop, the Hamiltonian of the 3q system in the rest frame has been derived. 
The methodology of the approach has been reviewed recently ^ and so will be sketched 
here only briefly. The Y-shaped baryon wave function has the form: 

By{xi, X2, X3, X) = eapyq'^ixi, X)q^{x2, X)q^{x3, X), (1) 

where a, /?, 7 are the color indices, q{xi,X) is the extended operator of the i^^ quark at 
a point Xi, and X = (0, X) is the equilibrium junction position. The wave function (fT|) 

^The preview of this calculation has been done in 0. 
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describes the only gauge-invariant configuration possible for baryons. The starting point 
of the approach is the Feynman-Schwinger representation for the gauge-invariant Green 
function of the three quarks propagating in the nonperturbative QCD vacuum 

3 °° 

G{x,y) = lljdsij Dz,exp{-K,){W)B, (2) 

1=1 Q 

where x = {xi,X2,X3}, y = {yi, ?/2, 1/3}, Zi = Zj(sj) are the quark trajectories with 
Zi{0) = Xi, Zi(T) = yi, while is the Fock-Schwinger proper time of the i^^ quark. 
Angular brackets mean averaging over background field. The quantities Ki are the kinetic 
energies of quarks, and all the dependence on the vacuum background field is contained 
in the generalized Wilson loop W 

W = ^_e,,keimnUi'urut, (3) 

with 

Uk = Pexpitg [ A^{x)dx''), k = 1,2,3. (4) 



Here P denotes the path-ordered product along the path Ci in Fig. 1 where the contours 
run over the classical trajectories of static quarks. In this figure three quark lines start at 
junction X at time zero, run in the time direction from to T with the spatial position 
of quarks fixed and join again in the junction Y at time T. There are three planes that 
are bounded, respectively, by one quark line, two lines connecting junction and quark at 
t = and t = T, and the connection line of two junctions. Under the minimal area law 
assumption, the Wilson loop configuration takes the form 

{W)b oc exp(-a(Si + ^2 + ^3)), (5) 

where Si are the minimal areas inside the contours formed by quarks and the string 
junction trajectories and a is the QCD string tension. 

In Eq. (0) the role of the time parameter along the trajectory of each quark is played by 
the Fock-Schwinger proper time Sj. The proper and real times for each quark are related 
via a new quantity that eventually plays the role of the dynamical quark mass. The final 
result is the derivation of the Effective Hamiltonian (EH), see Eq. (jH)) below. 

In contrast to the standard approach of the constituent quark model, the dynamical 
masses rrii are no longer free parameters. They are expressed in terms of the running 
masses ■mf\Q'^) defined at the appropriate hadronic scale of from the condition of 
the minimum of the baryon mass as a function of m^. 

Technically, this has been done using the einbein (auxiliary fields) approach, which is 
proven to be rather accurate in various calculations for relativistic systems. Einbeins are 
treated as c-number variational parameters: the eigenvalues of the EH are minimized 
with respect to einbeins to obtain the physical spectrum. Such procedure, first suggested 
in 0, [Hj, provides the reasonable accuracy for the meson ground states 

The method was already applied to study baryon Regge trajectories ^ and very recently 
for computation of magnetic moments of light baryons ^U] • The essential point adopted 
in j3] and continued in this paper is that it is very reasonable that the same method 
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should also give a unified description for both light and heavy baryons including baryons 
with two heavy quarks. As in IQ,^ we take as the universal QCD parameter the string 
tension a. We also include the perturbative Coulomb interaction with the frozen strong 
coupling constant a^. 

From experimental point of view, a detailed discussion of the excited QQ'q states is 
probably premature. Therefore we consider the ground-state baryons without radial and 
orbital excitations, in which case tensor and spin-orbit forces do not contribute pertur- 
batively. Then only the spin-spin interaction survives in the perturbative approximation. 
In what follows we disregard the hyper-fine splitting among the baryon masses, then the 
EH has the following form: 

Here Hq is the nonrelativistic kinetic energy operator and V is the sum of the perturbative 
one-gluon-exchange potential Vc'- 
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i<3 ■> 



where Vij are the distances between quarks, and the string potential String- The baryon 
mass is given by formula 

Mb = mm{H) + C, (7) 



where C is the quark self-energy correction calculated in [TI] : 



with 7] = 1 for g-quark^, t] = 0.88 for s-quark, t] = 0.234 for c-quark, and t] = 0.052 for 
6-quark. 

The string potential calculated in [7j as the static energy of the three heavy quarks was 
shown to be consistent with that given by a minimum length configuration of the strings 
meeting in a Y-shaped configuration at a junction X: 

Ktrmg(^^l, rg) = Cr/min, (9) 

where /min is the sum of the three distances between quarks and the string junction point 
X. The Y-shaped configuration was suggested long ago [T2j, and since then was used 
repeatedly in many dynamical calculations [13]. 



3 Solving the three-quark equation 
3.1 Jacobi coordinates 

The baryon wave function depends on the three-body Jacobi coordinates 

P., = Min-r^ (10) 

V A'' 



^Here and throughout the paper q denotes a hght quark u or d 
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^ !!^,r^±^_^\ (11) 

]j fj, \ mi + nij J 

{i,j,k cyclic), where /lij and /lij^k are the appropriate reduced masses 

_ rnjmj _ {■mi + mj)mk , . 

mj + rrij rrii + rrij + irik 

and iJ, is an arbitrary parameter with the dimension of mass which drops off in the final 
expressions. The coordinate p^j is proportional to the separation of quarks i and j and 
coordinate Xij is proportional to the separation of quarks i and j, and quark k. There 
are three equivalent ways of introducing the Jacobi coordinates, which are related to each 
other by linear transformations with the coefficients depending on quark masses, with 
Jacobian equal to unity. In what follows we omit indices 

In terms of the Jacobi coordinates the kinetic energy operator Hq is written as 



where R is the six-dimensional hyperradius, 

R^^ + (14) 

Q denotes five residuary angular coordinates, and K'^{D,) is an angular operator whose 
eigenfunctions (the hyperspherical harmonics) are 

K\Q)Y[K] = -K{K + 4)Y[K], (15) 

with K being the grand orbital momentum. In terms of Y^k] the wave function ip{p, A) 
can be written in a symbolical shorthand as 

K 

In the hypcrradial approximation, which we shall use below, K = and ip = ip{R). 
Obviously, this wave function is completely symmetric under quark permutations. Note 
that the centrifugal potential in the Schrodinger equation for the reduced radial function 
X{R) = R^/'^iPk{R) with a given K 

{K + 2f - 1/4 

is not zero even for K = 0. 

The Coulomb potential is directly expressed in terms of Jacobi coordinates: 

2as ^ [Jiij 1 

i<j 



while the corresponding expression for the string potential is more complicated. We will 
construct it in the next section. 
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3.2 String junction point 



Consider the definition of tlie minimal lengtfi string Y-sliaped configuration. Let (pijk be 
the angle between the line from quark i to quark j and that from quark j to quark k. 
One should distinguish two cases. If all ipijk are smaller than 120°, and the equilibrium 
junction position X coincides with the so-called Torrichelli point of the triangle in which 
vertices three quarks are situated. If ipijk is equal to or greater than 120°, the lowest 
energy configuration has the junction at the position of quark j. 

To find the Torrichelli point, consider a scalar function C{r) of a point inside a triangle 
AABC, defined as a sum of distances between this point and the triangle vertices: 

£(r) = |r - r^l + |r - r^l + |r - rc\- (17) 

djC' 

The position of the minimum of the function C is calculated from the condition — — = 0, 

dr 

i.e.: 

r - va r - rs r - rc „ 

1 r + 1 r + 1 r = -riA - - nc = 0, (18) 

\r — rA\ \r — rB\ \r — rc\ 

where nA,B,c are the unit vectors from the Torrichelli point directed to the vertices of 
the triangle. It follows from equation (fTH|) that this condition can be realized only in 
the case when angles between vectors n^i, n^, and ric are equal to 120°. If ipijk are all 
smaller than 120°, this point exists and is the unique one. From this point all sides of 
the triangle are seen at angle of 120°. 

The geometrical construction of the Torrichelli point is presented on Fig. 2. One should 
plot three equilateral triangles AAFB, ABDC, ACEA on the sides of the initial triangle 
AABC. It is lightly to prove the following statements: 

• The straight lines AD, BE, CF are crossed at a unique point T; 

• AD = BE = CF = AT + BT + CT; 

• ZATF = ZBTF = ZBTD = ZCTD = ZCTE = ZATE = -. 

3 

Now one can easy obtain an expression for a radius-vector of the Torrichelli point in 
terms of the lengths U of the segments between this point and the i^^ quark (segments 
AT, BT, CT on Fig. 2), and the quark positions Vi (points A, B, C on Fig. 2) jHj : 

^2^3 + ^1^3 + hh 

An equivalent expression for X in terms of the center-of-mass position -Rem, and vectors 
p and A is ^3] 

X = i?cm + «P + P\ (20) 



with 

a 




mi 



1 4t+ (3-^2) cot X 



rrii + rrij y 3/ijj 2A sin x 1 + t'^ 
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where 

2Asinx + W P 

^ ^ V 



2ACOSX+ X — — P 

p,ij rrii + rrij 



and X is the angle between p and A, < x < ti" ■ It can be easily seen that dependence 
on rrii in Eq. (jSOI) is apparent and X does not depend on quark masses just as it should 
be. 

Eq. is not very useful for our purposes. What we really need is the explicit expression 
for /^jjj in terms of the Jacobi coordinates [12]. Introducing the variable 9 = arctan(p/A), 
< < 7r/2 one obtains for the case ipijk < 120°: 

^min = "^os^ 9 X 

/ (m? — mn) tan^ 9 ( mo — mi ^ \ tan 9 1 \ , , 

X r ' . y 2Y+ ' ' cosx + v^sinx + , 21 

\mim2{mf - 7712) \m2 + mi J m P-i2,3 J 

where = mim2m,^/{mi + m2 + ^3). If mi = m2 = ms, this expression coincides with 
that derived in [Zj. 

If (^ijk > 120°, /min = Tij + Tjfc, whcre 



ri2 =Rsm9 ' ^ 



P'12,3 



„ I P' / m^ 2m 

ri3 =it cos 9 ■ . ■ \ / — 5- tan fc* H tan 6* cos x + 1 ( 22) 

AH2,3 V mf mi ^ ^ 



/ /"^ / 9 /I 2m 
?"23 =-R cos 9 ■ . ■ \ / — 5- tan t' tan 9 cos x + 1 

V /"12,3 V "^2 "^2 

The boundaries corresponding to the condition (pijk = 120° in the (x, 9) plane are: 

Oi(2) ix) = arctan(mi(2) (=F cos x - sin x/ V^)/m), ^^^^ 
93{x) = arctan(m2(/(x) + + 4x)/2m), 

where 

/(X) = (1 - x) cos X + (1 + x) sin x/ V3, (24) 

and X = mi/m2. These boundaries are shown in Fig. 3 for the case of equal quark 
masses. 

For simplicity, the string junction point is often chosen as coinciding with the center-of- 
mass coordinate. In this case 

String = O" ^ ] ■ y/jHHjji ' (25) 

{i,j,k cyclic). Accuracy of this approximation that greatly simplifies the calculations 
was discussed in [Zj, [TH]. We shall comment on this point later on. 
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3.3 Hyperradial approximation 



Introducing the variable x = y/JtR and averaging the interaction U = Vc + V^tring over 
the six-dimensional sphere Qq, one obtains the Schrodinger equation for x{^)'- 



+ 2 



' a ^ 15 
+ - - 6x - — 

X QX-^ 



X{x) = 0. (26) 



Because the wave function if) must be finite at the origin x{x) ~ 0{x^^'^) as x ^ 0. As 
X ^ oo one can neglect the Coulomb-like and centrifugal terms, and Eq. (j^Bj) becomes 

^^-ZXW = 0, z = {2bfl'x. (27) 
This is the familiar Airy equation whose solution Ai(2;) behaves at infinity as 

ki{z) r^-Ti-^/^z~^/^e^^{--z''/^), Rez> 0. (28) 
2 3 

In Eq. (jSni) Eq is the ground-state eigenvalue and 

a = Ry/Jl ■ j Vc{ri, rs, rg) ■ cifig = ' ^ ' ' 



h = 



^-=1 ■ j Ktring(T'i, r2, rs) ■ dile. 



<i / (29) 



The expression for the coefficient h can not be obtained analytically except for the equal 
quark mass system mi = m2 = = m, in which case the straightforward calculation 
yields: 

4 a f ^ I2V2 3^3 l\ a 

b=—-^-\ -VS + + ^— arccos - ^ 1.58 ■ 30 

15 Vm \ TT TT 5 / \/m 



On the contrary, in the approximation ()25|) the coefficient b can be easily found analyti- 
cally for the case of arbitrary quark masses: 



b = a-^- y^^]. (31) 
157r \ rrifc ' 



Let us explain the numerical coefficients in (jSHI), in more details. To this end we 
introduce the angle 6 as in ()2ip . such that 

p = Rsme, X = Rcose, 0<9<^, (32) 

and write the volume element d^pd^X as 

rf3p _ {A-nfp'^X^dpdX = {A-nfR^ sin^ 6 cos^ OdRdO. (33) 
The volume of the six-dimensional sphere is 

TT 

2 

f]g = (47r)2 /" sin^ e cos^ OdO = {Aii f ■ — = tt^. (34) 
J 16 
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Then averaging the Coulomb and string terms yields 



2 



< i >= 4 ■ - ■ (4vr)' ■ [ sin^cos^^ d9 = — ■ (35) 
p 71'^ R J Svr it 







2 



< A >= ■ i? ■ (47r)2 • / sin^ 9 cos^ 6 dO = ^ ■ R. (36) 





Combining together expressions §^ lead to (1^ . 



3.4 Analytic results for light baryons 

We can eliminate all dimensional parameters from the equation (j2EI) by a substitution 
y = b^/^x, which leads us to the equation: 

where 

The eigenvalue of the Eq. ()37p can be found using the ordinary perturbation theory, 
considering the Coulomb term f — ) as a small perturbation. This approximation works 

V yj 

well for a nucleon containing three light quarks with the running mass equal to zero. In 
this case there is only one dynamical quark mass m. So, the task is greatly simplified 
and one can obtain analytic expressions for m and Mb via two parameters: a and a^, as 
expansions in powers of as- 



Omitting the intermediate details outlined in Appendix we quote here the result with 
accuracy up to a^: 

m = 0.959 ■ • (1 + 0.270a, + O.llTa^ + ...), (38) 
Mb = 5.751 ■ ■ (1 - 0.270a, - 0.039a^ + ...) + C, (39) 

or 

Ms = 6m ■ (1 - 0.540a, - 0.083a^ + ...) + C, (40) 

where the constant C has been defined in Eq. (jH}. As it follows from Eq. ()40|) . the 
Coulomb-like correction to Mb comprises approximately 20%. 



3.5 Quasi-classical solution 

For the purpose of illustration, the problem is first solved quasi-classically rather than 
using quantum mechanics. This approach is based on the well known fact that interplay 
between the centrifugal term and the confining potential produces an effective potential 
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minimum specific for tlie tliree-body problem. Tfie numerical solution of (j26|) for the 
ground-state eigenenergy may be reproduced on a percent level of accuracy by using the 
parabolic approximation ^\ for the effective potential 

U{x) = — + bx + —. 

This approximation provides an analytical expression for the eigenenergy. The potential 
U{x) has the minimum at a point x = xq, which is defined by the condition U'{xq) = 0, 
i.e.: 

bxl + axo - 15/4 = . (41) 
Expanding U (x) in the vicinity of the minimum one obtains 

U{x) ^ U{xo) + ^U"{xo){x - xo)\ 

i.e., the potential of the harmonic oscillator with the frequency lo = ^JU"{xq). Therefore 
the ground-state energy eigenvalue is 

Eo^f/(xo) + ^^. (42) 



3.6 Variational solution 



Another rather accurate method of solving Eq. ()26|) is the minimization of the baryon 
energy using a simple variational Ansatz 

X(x) ~ x^/2e-^'^^ (43) 

where p is the variational parameter. Then using the three-quark Hamiltonian one can 

get an approximate expression for the ground-state energy: Eq min£'o(p), where 

p 

„ / X , , , V o 3 /TT , 15 pK 1 

Eo{p) = {x\H\x) = Sp'-a--- ^--p + b--- (44) 
The accuracy of the variational solution will be illustrate in Sect. 4. 



3.7 Analytical results for (Qud) baryons 

For the heavy quarks (Q=c and b) the variation in the values of their masses tuq is 
marginal. This is illustrated by the simple analytical results for (Qud) baryons jT5j . 
These results were obtained from the approximate solution of equation 

(9Eo(mi,m2,m3,p) ^ ^ ^^^^ 
dp 

where Eq is given by Eq. ()44j) in the form of expansion in the small parameters 

and as, (46) 



(0) 
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where rnq' is the heavy-quark running mass. 
Omitting the intermediate steps one obtains: 

Eo = Sv^^y [l + A-i~^-B-as + ..^, (47) 

^9 = ^[-) {^-A-i + B-as + ...), (48) 

mQ = mg)(l + O(e^a^,a.0 + •••)' (49) 
where for the Gaussian variational Ansatz 



Note that the corrections of the first order in ^ and are absent in the expression PU]) 
for mq. Accuracy of this approximation is illustrated in Table 1 of [Tf 



4 Baryon masses 

4.1 Quark dynamical masses 

We first calculate the dynamical masses irii retaining only the string potential in the 
effective Hamiltonian (jHl). This procedure is in agreement with the strategy adopted in 
[in] . The masses are then obtained from the equation: 

= 0> (51) 

orrii 

where 

5^ ( ^^ + ^ ) +^0(^1,^2, m3;a, = 0). (52) 

Then we add the one-gluon-exchange Coulomb potential and solve Eq. to obtain 

the ground-state eigenvalues i?o('^i; "^2, ""^s; «s) for a given a^- The physical mass Mb 
of a baryon is 

Mb = 5^ I + Y 1 + + ^^^^ 




i=l 



We use the values of parameters a = 0.15 GeV^ (this value has been found in a re- 
cent lattice study [IH]) , = 0.39, mf^ = 0.009 GeV (g = u,d), mi°^ = 0.17 GeV, 
m^c^ = 1.4 GeV, and m[°'' = 4.8 GeV. The results for various baryons, obtained using 
various approximations, are given in Tables 1-3. Table 1 contains the results obtained us- 
ing the variational solution of Eq. using the approximation for the three-quark 
potential . In Table 2 the results are shown obtained using the same approximation. 
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X = Rem, but exact numerical solution of Eq. (|26|) . Table 3 contains the results ob- 
tained by the numerical integration of with the genuine three-quark potential in the 
form (j2H)j (j22I)' Comparing the results of Tables 1 and 2 we observe a good accuracy 
of the variational solution: the difference between variational and exact results for Mb 
does not exceed 10 — 15 MeV for all baryons from lightest to doubly heavy ones. The 
approximation X = Rem leads to a ~ 5% increase of the coefficient b in (j26|l . As a 
consequence, this approximation increases the baryon masses by ~ 70 MeV (compare 
the results of Tables 2 and 3. 

Note that there is no good theoretical reason why quark masses rrii need to be the same 
in different baryons. Inspection of Table 1 shows that the masses of the light quarks 
{u, d, or s) are increased by ~ 100 MeV when going from light to heavy baryons. The 
dynamical masses of light quarks rrtu^^s ~ ~ 400 — 500 MeV qualitatively agree with 
the results of [20j obtained from the analysis of the heavy-light ground-state mesons. 

While studying Table 3 is sufficient to have an appreciation of the accuracy of our predic- 
tions, few comments should be added. We expect an accuracy of the baryon predictions 
to be ~ 5 — 10% that is partly due to the approximations employed in the derivations 
of the Effective Hamiltonian itself and partly due to the error associated with the 
variational nature of hyperspherical approximation. From this point of view the overall 

agreement with data is quite satisfactory. E.g. we get -(A^ + A)thcory = 1144 MeV vs. 
-(A^ + A)cxp = 1085 MeV ( a 5% increase in as would correctly give the (A^, A) center 

of gravity), ^(A + S + 2S*) = 1242 MeV vs. experimental value of 1267 MeV. We 

also find Stheory = 1336 MeV (without hyperfine splitting) vs. Sexp = 1315 MeV and 
■^theory = 2542 McV VS. Sg^p = 2584 MeV. On the other hand, our study shows some 
difficulties in reproducing the il-hyperon mass. 



4.2 Doubly heavy baryons 

In Table 4 we compare the spin-averaged masses (computed without the spin-spin term) 
of the lowest doubly heavy baryons to the predictions of other models j^I], [21] as well 
as variational calculations of j3] for which the center-of-gravity of non-strange baryons 
and hyperons was essential a free parameter. Most of recent predictions were obtained 
in a light quark-heavy diquark model ^2], in which case the spin-averaged values 
are M = |(Mi/2 + 2M3/2). Note that the wave function calculated in the hyperspheri- 
cal approximation shows the marginal diquark clustering in the doubly heavy baryons. 
This is principally kinematic effect related to the fact that in this approximation the 
difference between the various mean values f^j in a baryon is due to the factor ^JlJJi^j 
which varies between for rrij = nij and yjxjrrii for vrii <C vrij. In general, in 

spite of the completely different physical picture, we find a reasonable agreement within 
100 MeV between different predictions for the ground-state masses of the doubly heavy 
baryons. Our prediction for M^cu is 3.66 GeV with the perturbative hyperfine splitting 
S*^^ — Seen ~ 40 MeV. This splitting has been calculated using the Fermi-Breit spin-spin 
interaction [21]. The change of o to 0.17 GeV^ increases the mass of Sec by ~ 30 MeV. 
Note that the mass of Sec is rather sensitive to the value of the running c-quark mass 
mi"-*, see Fig. 4. 
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5 Conclusions 



In this paper, we have outhned a novel approach to baryon spectroscopy which is based on 
a single framework of the Effective Hamiltonian that is consistent with QCD. This model 
uses the stringlike picture of confinement and perturbative one-gluon-exchange potential. 
The main advantage of this work is demonstration of the fact that it is possible to describe 
all the baryons in terms of the only two parameters inherent to QCD, namely a and as- 
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Appendix 

The approximate calculation of the mass of the three-light-quark system in 

the hyperspherical formalism 



Let us consider the case when each quark has a zero current mass and the constituent 
mass m, which is the same for all three quarks. Then the mass of the system is the 
minimum of the function ?i(m): M. = min?i(m), where^ 

m 

tin 

H{m) = ^ + E, (54) 
and E is the energy level defined from the ordinary Shrodinger equation (see Eq. 

(fx 0/^7-1 PV^ 7 15 



+ 2[E + ^ ^^ = 0, (55) 



where 



13 = — ^ 2.401 ■ a„ (56) 

OTT 

7 = ^^a^ 1.663 -a. (57) 
157r 



Let us replace a variable x by a dimensionless variable y. y = { — x. Then: 

\m J 



^For simplicity we omit here and below the corrections C. 
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where £ and 5 are the dimensionless parameters: 

/ \ 1/3 / 2\ 1/3 

We will calculate the eigenvalues of the equation using the second-order perturbation 
theory, considering the Coulomb term ( — — ) in the potential as a small perturbation. 

V yJ 

The unperturbed equation is: 

It contains no physical parameters, so its solution is a pure mathematical task. Let us 
denote the eigenvalues of the equation as {Aj}: 

< Ao < Ai < . . . , 

and the corresponding normalized eigenfunctions as {fi{y)}- In what follows we will use 
the notations: 

+ 00 

e = j ^fo{y)dy, e>o, 



n 

+ 00 

J -My)fo{y)dy 

^ = ^ a;^^ ' 

i=l * ^ 

The approximate numerical values of these parameters are: 

Ao ^ 3.030, (60) 
^ ^ 0.553, (61) 
ri ^ 0.028. (62) 

The ground level of the equation can be approximately calculated as follows: 

8 ^ Xo-6^-6^r]. (63) 

The small parameter here is the ratio -p. To estimate this quantity one can solve the 

Ao 

task in the zero approximation without the Coulomb term. So, 

£ ~ Ao, 

E ^ Xo(-) , 

3 . fl'^'^' 



n ^ -m + Xo . 

2 \m 

Ti has a single minimum, defined by an equation dTi/dm = 0, i.e., 

I - IW^'m-'/' = 0. 
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Thus m = I V^' finally: 



Ao Ao V 7 y 



0.360a, 



For as = 0.4 one has — ~ 0.14 <^ 1. This verifies the correctness of using the perturba- 

Ao 

tion theory in this task. 

Now, using Eq. (jU^ . we can calculate the constituent mass m and the mass of the state 
Ai. For the energy level we have: 

/ 2\l/3 / 2\l/3 

^=(-) ^=(-) >^0-{imY/''p^-mP^ri. 

It is convenient to use the positive dimensionless parameter q = m^/^7~^/^, so that 
m = ^/^q^■ Then 

and substituting it to the Eq. we find: 

n = Vl (^oq-' - +(\- q' 



dl~L 

Ti has a single minimum, defined from the condition: — — 

oq 



= 0, i.e., 

q=go 



(64) 



After calculating qo, one can find the constituent mass m = ^J^qQ and the mass of the 
system M.: 

M = niqo) = M _ p^A . 

3 V 5o / 

The equation can be lightly solved: 



'' = ^1 9^W~V • ^''^ 

On expanding the right-hand side of this equation taking /3 as a small parameter 

we easily get formulas for the constituent mass: 
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and for the state mass: 



48Ao 6 



(67) 



Substituting in the equations and ^7\i numerical values according to formulas 
(jKTjl . p]|l . (jnH), dnSl), we obtain the results quoted in Eqs. (001 

For a = 0.15 GeV^ a, = 0.4 one has m ^ 0.418 GeV, M ^ 1.973 GeV. 
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Tables 



Table 1. Summary of variational calculations of Eq. (j^Uj) for the various baryon states 
in the approximation X = Rem- 



Baryon 




7712 


ms 




Mb 


(qqq) 


0.372 


0.372 


0.372 


1.433 


1.221 


(qqs) 


0.377 


0.377 


0.415 


1.404 


1.314 


(qss) 


0.381 


0.420 


0.420 


1.377 


1.405 


(sss) 


0.424 


0.424 


0.424 


1.350 


1.493 


(qqc) 


0.424 


0.424 


1.464 


1.178 


2.538 


(qsc) 


0.427 


0.465 


1.467 


1.153 


2.613 


(ssc) 


0.468 


0.468 


1.469 


1.129 


2.686 


(qqb) 


0.446 


0.446 


4.819 


1.093 


5.909 


(qsb) 


0.448 


0.487 


4.820 


1.067 


5.978 


(ssb) 


0.490 


0.490 


4.821 


1.042 


6.046 


(qcc) 


0.459 


1.498 


1.498 


0.914 


3.712 


(sec) 


0.499 


1.499 


1.499 


0.890 


3.777 


(qcb) 


0.477 


1.524 


4.834 


0.793 


7.021 


(scb) 


0.517 


1.525 


4.834 


0.770 


7.083 


(qbb) 


0.495 


4.854 


4.854 


0.606 


10.260 


(sbb) 


0.534 


4.855 


4.855 


0.583 


10.318 



Table 2. The results obtained from the exact numerical solution of Eq. (j^Bj) in the 
approximation X = Rem for the various baryon states . 



Baryon 


mi 


m2 


ms 




Mb 


(qqq) 


0.372 


0.372 


0.372 


1.427 


1.212 


(qqs) 


0.376 


0.376 


0.415 


1.398 


1.306 


(qss) 


0.381 


0.419 


0.419 


1.370 


1.397 


(sss) 


0.423 


0.423 


0.423 


1.344 


1.485 


(qqc) 


0.424 


0.424 


1.464 


1.170 


2.530 


(qsc) 


0.426 


0.465 


1.466 


1.146 


2.604 


(ssc) 


0.467 


0.467 


1.468 


1.122 


2.677 


(qqb) 


0.445 


0.445 


4.820 


1.085 


5.900 


(qsb) 


0.448 


0.487 


4.820 


1.059 


5.970 


(ssb) 


0.488 


0.488 


4.820 


1.035 


6.037 


(qcc) 


0.458 


1.497 


1.497 


0.905 


3.702 


(sec) 


0.497 


1.498 


1.498 


0.882 


3.767 


(qcb) 


0.475 


1.523 


4.833 


0.784 


7.010 


(scb) 


0.515 


1.523 


4.837 


0.760 


7.072 


(qbb) 


0.490 


4.850 


4.850 


0.596 


10.245 


(sbb) 


0.530 


4.856 


4.856 


0.571 


10.303 
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Table 3. The same as in Table 2 but for an exact treatment of the string-junction point. 



Baryon 


mi 


1712 


ms 


Eq 


Mb 


(qqq) 


0.362 


0.362 


0.362 


1.392 


1.144 


(qqs) 


0.367 


0.367 


0.407 


1.362 


1.242 


(qss) 


0.371 


0.411 


0.411 


1.335 


1.336 


(sss) 


0.415 


0.415 


0.415 


1.307 


1.426 


(qqc) 


0.406 


0.406 


1.470 


1.142 


2.464 


(qsc) 


0.409 


0.448 


1.471 


1.116 


2.542 


{ssc) 


0.452 


0.452 


1.473 


1.090 


2.621 


(qqb) 


0.425 


0.425 


4.825 


1.054 


5.823 


(qsb) 


0.429 


0.469 


4.826 


1.026 


5.903 


(ssb) 


0.471 


0.471 


4.826 


1.000 


5.975 


(qcc) 


0.444 


1.494 


1.494 


0.876 


3.659 


(sec) 


0.485 


1.496 


1.496 


0.851 


3.726 


(qcb) 


0.465 


1.512 


4.836 


0.753 


6.969 


(scb) 


0.505 


1.514 


4.837 


0.729 


7.032 


(qbb) 


0.488 


4.847 


4.847 


0.567 


10.214 


(sbb) 


0.526 


4.851 


4.851 


0.544 


10.273 



Table 4. Comparison of various predictions for ground-state masses (in units of GeV) 
of doubly heavy baryons. 



Baryon 


Ref. dS] 


Ref. i 


Ref. |2B 


Ref. [21 


Ref. 123 


' — 'CC 


3.66 


3.69 


3.57 


3.69 


3.70 




3.73 


3.86 


3.66 


3.84 


3.80 




6.97 


6.96 


6.87 


6.96 


6.99 


^cb 


7.03 


7.13 


6.96 


7.15 


7.07 


^bb 


10.21 


10.16 


10.12 


10.23 


10.24 


^bb 


10.27 


10.34 


10.19 


10.38 


10.34 
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Figure captions 



Figure 1. Three-lobes Wilson loop. 

Figure 2. The geometrical construction of the Torrichelli point T for an arbitrary 
triangle AABC. The triangles AAFB, ABDC, ACEA are equilateral. 

Figure 3. The four regions in the (x, 9) plane corresponding to (pijk > 120° and 
fijk < 120° for the case of equal quark masses. 

Figure 4. Mass of Sec as a function of the running c-quark mass for a = 0.15 GeV^ 
and a = 0.17 GeV^. The masses are given in GeV. Bold points refer to the case 
m^c^ = 1.4 GeV. 
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B 



Fig.2. 
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Fig.3. 
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